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Measurement of the tt Production Cross Section in pp Collisions at y/s = 1.96 TeV 
using Lepton + Jets Events with Lifetime ^-tagging 
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We present a measurement of the top quark pair (tt) production cross section (cr t j) in pp colli- 
sions at yfs — 1.96 TeV using 230 pb _1 of data collected by the D0 experiment at the Fermilab 
Tevatron Collider. We select events with one charged lepton (electron or muon), missing transverse 
energy, and jets in the final state. We employ lifetime-based fe-jet identification techniques to fur- 
ther enhance the tt purity of the selected sample. For a top quark mass of 175 GeV, we measure 
a t i = 8. 61^5 (stat. + syst.) ± 0.6 (lumi.) pb, in agreement with the standard model expectation. 

PACS numbers: 13.85.Lg, 13.85.Qk, 14.65.Ha 

The top quark was discovered at the Fermilab Tevatron sions at ^fs = 1.8 TeV based on about 50 pb -1 of data per 
Collider by the CDF and D0 collaborations [lj in pp colli- experiment. The increased statistics and higher collision 
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energy of yfs = 1.96 TeV of Tevatron Run II allow more 
precise measurements of top quark properties, including 
its production and decay characteristics. Theoretical cal- 
culations performed within the framework of the stan- 
dard model (SM) predict the tt production cross section 
(a t t) with an uncertainty of less than 15% 0. A signifi- 
cant deviation from this prediction would signal the pres- 
ence of physics beyond the SM, such as tt resonant pro- 
duction 0] . The CDF and D0 collaborations have previ- 
ously reported measurements of o t t at \fs = 1.8 TeV |4T 
Recent measurements at y/s = 1.96 TeV by the CDF 5j 
and D0 collaborations agree with the SM prediction 
within their experimental uncertainties. 

In the SM, the top quark decays to a W boson and b 
quark with a branching ratio of w 100%. The lepton+jets 
final state results from the leptonic decay of one of the W 
bosons and the hadronic decay of the other. The event 
signature is one lepton with high transverse momentum, 
large transverse energy imbalance ($r) due to the unde- 
tected neutrino, and four jets, two of which result from 
hadronization of the b quarks. 

In this Letter, we report the measurement of a t i in the 
lepton (electron or muon) plus jets channel using &-jet 
identification (^-tagging) techniques exploiting the long 
lifetime of B hadrons. The data were collected by the D0 
experiment from August 2002 through March 2004, and 
correspond to an integrated luminosity of 226 ± 15 pb _1 
(229 ± 15 pb _1 ) in the electron (muon) sample. 

The D0 detector includes a tracking system, calorime- 
ters, and a muon spectrometer 0. The tracking system 
consists of a silicon microstrip tracker (SMT) and a cen- 
tral fiber tracker (CFT), both located inside a 2 T super- 
conducting solenoid. The tracker design provides efficient 
charged particle measurements in the pseudorapidity re- 
gion | ry | < 3 8]. The SMT strip pitch of 50-80 allows 
a precise reconstruction of the primary interaction ver- 
tex (PV) and an accurate determination of the impact 
parameter of a track relative to the P V 9] , which are the 
key components of the lifetime-based 6-jet tagging algo- 
rithms. The PV is required to be within the SMT fiducial 
volume and consist of at least 3 tracks. The calorime- 
ter consists of a central section (CC) covering \q\ < 1.1, 
and two end calorimeters (EC) extending the coverage to 
\r)\ w 4.2. The muon system surrounds the calorimeter 
and consists of three layers of tracking detectors and two 
layers of scintillators [10(. A 1.8 T iron toroidal mag- 
net is located outside the innermost layer of the muon 
detector. The luminosity is calculated from the rate for 
pp inelastic collisions detected using two hodoscopes of 
scintillation counters mounted close to the beam pipe on 
the front surfaces of the EC calorimeters. 

We select data samples in the electron and muon chan- 
nels by requiring an isolated electron with pr > 20 GeV 
and 1 77 1 < 1.1, or an isolated muon with px > 20 GeV 
and \r)\ < 2.0. More details on the lepton identification 
as well as trigger requirements are reported elsewhere [6J . 



In both channels, we require $t to exceed 20 GeV and 
not be collinear with the lepton direction in the trans- 
verse plane. These W boson candidate events must be 
accompanied by one or more jets with pt > 15 GeV 
and rapidity \y\ < 2.5 0- Jets are defined using a cone 
algorithm with radius A1Z = 0.5 [ll|. We classify the 
selected events according to their jet multiplicity. Events 
with 3 or > 4 jets are expected to be enriched in tt sig- 
nal, whereas events with only 1 or 2 jets are expected 
to be dominated by background. We use the former to 
estimate a t t, and the latter to verify the background nor- 
malization procedure. 

The main background in this analysis is the production 
of W bosons in association with jets (W+jets), with the 
W boson decaying leptonically. In most cases, the jets 
accompanying the W boson originate from light (u, d, 
s) quarks and gluons (W+light jets). Depending on the 
jet multiplicity, between 2% and 14% of W^+jets events 
contain heavy flavor jets resulting from gluon splitting 
into bb or cc (Wbb or Wcc, respectively), while in about 
5% of events, a single c quark is present in the final state 
as a result of the W boson radiated from an s quark 
from the proton's or antiproton's sea (Wc). A sizeable 
background arises from strong production of two or more 
jets ("multijets"), with one of the jets misidentified as 
a lepton and accompanied by large J5t resulting from 
mismeasurements of jet energies. Significantly smaller 
contributions to the background arise from single top, 
Z+jets, and weak diboson (WW, WZ and ZZ) produc- 
tion. Only a small fraction of the background events 
contain b or c-quark jets in the final state. As a con- 
sequence, the signal-to-background ratio is significantly 
enhanced when at least one jet is identified as a 6-quark 
jet. 

We use a secondary vertex tagging (SVT) algorithm 
to identify fe-quark jets. Secondary vertices are recon- 
structed from two or more tracks satisfying the following 
requirements: pt > 1 GeV, > 1 hits in the SMT lay- 
ers and impact parameter significance d ca /<Jd ca > 3.5 |9(. 
Tracks identified as arising from Kg or A decays or from 
7 conversions are not considered. If the secondary ver- 
tex reconstructed within a jet has a decay length signif- 
icance L xy /(TL xy > 7 0| > the jet is tagged as a 6-quark 
jet. Events with exactly 1 (> 2) tagged jets are referred 
to as single-tag (double-tag) events. We treat single-tag 
and double-tag events separately because of their differ- 
ent signal-to-background ratios. 

Secondary vertices with L xy /aL xy < —7 appear due 
to finite resolution of their characteristics after recon- 
struction, and define the "negative tagging rate". The 
negative tagging rate is used to estimate the probability 
for misidentifying a light flavor (u, d, s quark or gluon) 
jet as a 6-quark jet (the "mis-tagging rate"). 

We estimate both the 6-tagging efficiency and the mis- 
tagging rate using jets with > 2 tracks satisfying less 
stringent requirements than those for SVT. In particular, 
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FIG. 1: Measured 6-tagging efficiency (circles) and mis- 
tagging rate (triangles), and estimated c-tagging efficiency 
(solid line) as a function of jet pr- 



the pt cut is reduced from 1 GeV to 0.5 GeV for all 
but the highest px track, and no cut on d ca /<Td ca of the 
tracks is made. These requirements have an efficiency 
per jet > 80% for px > 30 GeV and integrated over 
y. We measure the 6-tagging efficiency in a data sample 
of dijet events with enhanced heavy flavor content by 
requiring a jet with an associated muon at high transverse 
momentum relative to the jet axis. By comparing the 
SVT and muon-tagged jet samples, the tagging efficiency 
for semileptonic 6-quark decays ( "semileptonic 6-tagging 
efficiency") can be inferred. We make use of a Monte 
Carlo (MC) simulation to further correct the measured 
efficiency to the tagging efficiency for inclusive 6-quark 
decays. We estimate the c-tagging efficiency from the 
same simulation, corrected by a scale factor defined as the 
ratio of the semileptonic 6-tagging efficiency measured in 
data to that measured in the simulation. We estimate the 
mis-tagging rate from the negative tagging rate measured 
in dijet events, corrected for the contribution of heavy- 
flavor jets and the presence of long-lived particles in light- 
flavor jets. Figure ^ shows the 6-tagging efficiency, c- 
tagging efficiency and mis-tagging rate as a function of 
jet pt- 

We simulate tt production, and all background pro- 
cesses except multijets, using ALPGEN Il3l to generate 
the parton-level processes, and PYTHIA |l4| to provide 
fragmentation and to decay all unstable particles ex- 
cept B hadrons and r leptons, which are modeled via 
evtgen and TAUOLA 0, respectively. We process 
the generated events through the full GEANT-based 01 
D0 detector simulation and the same reconstruction pro- 
gram used to process the data. We apply small addi- 
tional smearing to the reconstructed objects to improve 
the agreement between the data and the simulation, and 



account for remaining discrepancies using correction fac- 
tors derived by comparing the efficiencies measured in 
Z — ► £ + £~ data to those obtained from the simulation. 
For all processes except the multijets background, we 
make use of the MC simulation to compute the total 
acceptance, applying the trigger, reconstruction and tag- 
ging efficiencies measured using data. The tagging proba- 
bility for a particular process depends on the flavor com- 
position of the jets in the final state as well as on the 
overall event kinematics. We estimate it by applying the 
tagging rates measured in data to each jet in the sim- 
ulation, taking into consideration its flavor, pt, and y. 
In the case of VF+jets events, we also use the simulation 
to estimate the fraction of the different VF+heavy flavor 
subprocesses. 

We compute the tt acceptance for events with a true 
electron or muon arising from a W — > Iv (t = e, fi, r) 
decay, corresponding to total branching fractions of 
17.106% and 17.036% 0, respectively, in the electron 
and muon channels. In the electron channel, the total ac- 
ceptance before tagging is estimated to be (10.8 ± 0.8)% 
and (14.2 ± 1.7)%, for events with 3 and those with > 4 
jets, respectively. The corresponding numbers for the 
muon channel arc (9.9 ± 1.0)% and (14.1 ± 1.9)%. The 
estimated single-tag efficiencies are (43.4 ± 1.2)% and 
(45.3 ± 1.0)% for events with 3 and those with > 4 jets, 
respectively. The corresponding double-tag efficiencies 
are (10.4 ± 1.0)% and (14.2 ± 1.3)%. 

We estimate the number of multijet events from the 
data for each jet multiplicity using the matrix method 
described in Ref. [||, separately for the samples before 
and after tagging. Smaller contributions from single top, 
Z+jets, and diboson production (collectively referred to 
as "other bkg") are estimated from the simulation, nor- 
malized to the next-to-leading order theoretical cross sec- 
tions [n| Ell- We also include under "other bkg" the 
contribution from tt with both W bosons decaying lep- 
tonically, assuming the same a t i as for tt — > £+jets. We 
determine the number of tagged VF+jets events as the 
product of the number of VF+jets events in data before 
tagging and the average tagging probability for VF+jets 
events (e.g. ~ 4% for single-tag and « 0.4% for double- 
tag events with > 4 jets). The number of W+jets events 
before tagging is computed as a difference between the 
number of selected events and the estimated contribu- 
tion from the rest of processes (multijets, tt, and "other 
bkg"). 

Tables HI and El summarize the sample composition for 
single-tag and double-tag events, respectively, assuming 
a t i — 7.0 pb. Figure |2 shows the observed and expected 
number of events for each jet multiplicity. We interpret 
the excess over the background expectation in the third 
and fourth jet multiplicity bins as the tt signal. The 
good agreement between observation and expectation in 
the first and second jet multiplicity bins validates the 
background estimation procedure. 



6 



c 

> 



n 
E 

3 



300 



200 



100 



(a) 



• D0data 

□ W+light 

O Wc,Wcc,Wbb 
■ Multijets 

□ Other bkg 

□ tf -> l+jets 



3 >4 
Jet multiplicity 



CO 

c 

5 30 



0) 

E 

3 



20 



10 



(b) 



• D0data 
□ Wbb 
I I Other bkg 
I I tt -» l+jets 



3 >4 
Jet multiplicity 



FIG. 2: Expected and observed number of (a) single-tag and (b) double-tag events. The hatched area represents the total 
uncertainty in the expectation. 



TABLE I: Summary of observed and expected numbers of 
events before tagging and with exactly one jet tagged. 





W+ljet 


IV+2jets 


VK+3jets 


W+> 4 jets 






Before 


tagging 




Observed 


14054 


5502 


1365 


367 


Multijets 


718±78 


516±43 


190±14 


66±6 


After tagging 


W+light 


36.8±4.0 


21.4±2.4 


7.2±0.9 


1.8±0.3 


Wc 


47.8±5.4 


24.2±2.7 


5.7±0.7 


0.8±0.1 


Wee 


12.2±3.4 


17.2±4.8 


6.6±1.9 


2.2±0.7 


Wbb 


33.9±8.7 


43.2±11.0 


15.1±3.9 


4.5±1.3 


Multijets 


14.9±1.9 


16.3±2.1 


8.3±1.5 


4.0±1.2 


Other bkg 


6.6±0.9 


20.6±2.1 


8.2±0.8 


2.2±0.4 


Total bkg 


152.4±14.8 142.9±16.0 


51.0±5.6 


15.6±1.7 


tt — * ^+jets 


0.4±0.1 


6.8±1.4 


24.4±1.7 


34.8±4.3 


Total expected 


152.8±14.8 149.7±16.2 


75.4±5.9 


50.4±4.8 


Observed 


137 


145 


79 


62 


TABLE II: Summary of observed and expected number of 


events with two 


or more jets tagged. 








W+2jets W- 


-3jets 


W+>4jets 


Wbb 


5.7±1.6 2.2±0.6 


0.7±0.2 


Other bkg 


3.7±0.4 2.0±0.3 


0.5±0.3 


Total bkg 


9.4±1.8 4.2±0.8 


1.2±0.3 


tt -> £+jets 


0.8±0.2 5.9±0.7 


10.9±1.9 


Total expected 


10.2+l.c 


10.1±1.2 


12.1±2.0 


Observed 


15 




9 


14 



We calculate a t i by maximizing a likelihood function 
including a Poisson term for each of the eight indepen- 
dent channels considered: 3 and > 4 jets, for single- 
and double-tag events in the electron and muon chan- 
nels. At each step in the maximization, the multijet 
background in these eight tagged samples, and the cor- 
responding samples before tagging, is constrained within 



errors to the amount determined by the matrix method. 
In addition, we include a Gaussian term for each of the 
systematic uncertainties considered, following the proce- 
dure described in Ref . |2l| . In this approach, each source 
of systematic uncertainty is allowed to affect the central 
value of the cross section during the maximization proce- 
dure, thus yielding a combined statistical and systematic 
uncertainty on <r t j. Assuming a top quark mass (m t ) of 
175 GeV, we measure 

o tl = S.6±i J (stat. + syst.) ± 0.6 (lumi.) pb, 

in good agreement with the SM prediction of 6.77 ± 0.42 
Pb 0. 

The contribution due to each individual source of sys- 
tematic uncertainty can be estimated by redoing the fit 
after fixing all but the corresponding Gaussian term and 
unfolding the statistical uncertainty from the resulting 
total uncertainty. The statistical uncertainty of pb 
is obtained from the fit where all Gaussian terms are 
fixed. As shown in Table IIIII 6-jet tagging efficiency, 
jet energy calibration, and background modeling are the 
leading sources of systematic uncertainty. In addition, a 
systematic uncertainty of 6.5% from the luminosity mea- 
surement [2^] has been assigned. In the top quark mass 



TABLE III: Summary of systematic uncertainties on a t t- 


Source 


Ao- t t (pb) 


6-tagging efficiency 


+0.6 -0.5 


Jet energy calibration 


+0.5 -0.4 


Background modeling 


±0.5 


Lepton selections 


+0.5 -0.4 


Jet identification 


+0.3 -0.2 


Multijet background 


+0.3 -0.2 


Mis-tagging rate 


±0.1 


Total 


+1.1 - 1.0 



7 



range of 160 GeV to 190 GeV, the measured cross section 
decreases (increases) by 0.06 pb per 1 GeV shift of m t 
above (below) 175 GeV. 

We used an alternative 6-tagging algorithm to cross 
check this result. This algorithm relics on counting tracks 
with significant impact parameter (CSIP) with respect to 
the PV: a jet is tagged if > 2 (> 3) associated tracks 
have d ca /od ca > 3 {d C a/(?d aa > 2). As compared to 
SVT, this algorithm has a slightly higher 6-tagging effi- 
ciency and about a factor of two higher mis-tagging rate. 
The measured cross section using the CSIP algorithm is 
a t t = 7.6l^4 (stat. + syst.) ± 0.5 (lumi.) pb, consistent 
with the SVT result once the existing overlap between 
both samples is taken into account. While we are cur- 
rently not combining these two results, the fact that dif- 
ferent 6-tagging techniques are only partially correlated 
will be exploited in future analyses to further increase 
the precision of this measurement. 

In summary, we have measured the tt production cross 
section in pp interactions at ^/s — 1.96 TeV in the lep- 
ton+jets channel using lifetime 6-tagging. Our measure- 
ment yields er t j = 8. 6+^5 (stat. + syst.) ± 0.6 (lumi.) pb, 
in a good agreement with the SM prediction. 
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